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I show how to protect adiabatic quantum computation (AQC) against decoherence and certain 
control errors, using a hybrid methodology involving dynamical decoupling, subsystem and stabilizer 
codes, and energy gaps. Corresponding error bounds are derived. As an example I show how to 
perform decoherence-protected AQC against local noise using at most two-body interactions. 



Adiabatic quantum computation (AQC), originally de- 
veloped to solve optimization problems offers a fas- 
cinating alternative to the standard circuit model [2| to 
which it is computationally equivalent [3|. The effects 
of decoherence on AQC were studied in several works 
0> S @| ■ Unlike the circuit model, for which an elabo- 
rate theory of fault tolerant QC exists along with a noise 
threshold for fault tolerance |7J, it is not yet known how 
to make AQC fault tolerant. Here I show how AQC can 
be can be protected against decoherence and certain con- 
trol errors. To do so I devise a hybrid method that involv- 
ing dynamical decou plin g (DD) [§] , subsystem S 1(J 11 1 
and stabilizer codes 12], and energy gaps [H,[lJ|- 



Viewed as a closed system, AQC proceeds via slow evo- 
lution on a timescale set by the system's minimal energy 
gap A from the ground state [H, 0| . In the presence of 
the system-bath interaction Hsb this gap can be signifi- 
cantly reduced because the interaction will cause energy 
level splittings, or an effective broadening of system en- 
ergy levels; when these levels overlap adiabaticity breaks 
down and so does AQC, even at zero temperature [5j. A 
bath at finite temperature presents another problem: in 
the universality proofs [3( the system energy gap scales 
as an inverse polynomial in the problem size, so that the 
temperature too must be lowered polynomially to prevent 
thermal excitations. All of the problems listed above are 
due to the presence of Hs b ■ Clearly, if Hsb can be effec- 
tively eliminated or reduced, this will enhance the fidelity 
of AQC. The main tool I shall use to this end is dynam- 
ical decoupling, which involves the application of strong 
and fast pulses. Perhaps surprisingly, this can be done 
without interfering with the slow adiabatic evolution. 

Distance measure and operator norm. — As a dis- 
tance measure between states I use the trace distance 
D[ Pl ,p 2 ] = HIpj - paid, where ||A||i = Tr|A|, \A\ = 
VaTa 0. When applied to pure states pi — 
I shall write -D[0i, V^]- The operator norm is ||A|| = 
sup|||^,)|| = i \\A\ip)\\ = max, A i; where A; E Spec(|A|). 

Closed-system adiabatic error. — Let s = t/T € [0, 1] 
be the dimensionless time, with T the final time. Let 
the system Hamiltonian that implements AQC, H at \{s), 
act on n qubits. In AQC the ground state |</> a d(s)) 
of iJad(s) at the final time s = 1 encodes the solu- 
tion to the computational problem [l[. The actual final 



state IV^l)) is the solution of the Schrodinger equation 
d\ip) /ds = — iTi? a d|^) {h — 1 units are used throughout). 
In AQC one is therefore interested in minimizing the er- 
ror <5 a d = D{ip(l), </> a d(l)]. Most of the known AQC algo- 
rithms interpolate between initial and final local Hamil- 
tonians, H and Hi, via H ad (s) = (1- f(s))H + f(s)H 1 , 
where /(0) = and /(l) = 1, and exhibit a final time 
that scales as a polynomial in the problem/system size 
n. Locality means that ||-ff a d|| ~ A O(n), where A is 
the energy scale. Thus \\d j H^/dsiW ~ A \d j f /ds j \0{n). 
Let {Ei(s)}i=o be the eigenvalues of i? a d(s), and let 
A = min, )S \Ei(s) — Eq(s)\ be the minimum gap from 
the instantaneous ground state energy E (s). Assume 
that A(n) ~ Aon~ 2 , where z > is the dynamical criti- 
cal exponent. Depending on the differentiability of i/ a d, 
and assuming that £f a d(0) = i? a d(l) = 0, one can prove 
different versions of the adiabatic theorem. For exam- 
ple, (i) [15]: if -ff a d(s) is twice differentiable on [0, 1] then 
provided T ~ r||i? a d|| 2 /A 3 the error can be made arbi- 
trarily small in the time dilation factor r > 1: <5 a d < r~ 2 . 
Or, (ii) [l6j |: if H a d(s) is infinitely differentiable on [0, 1] 
then provided T ~ rAf||7j a d||/A 2 , the error can be made 
exponentially small in the order N of an asymptotic ex- 
pansion: (5 a d < 7'^. In both cases 



T 



(1) 



where ( — 'iz + 2 for case (i) and £ = 2z + 1 for case 
(ii), and I omitted \d 3 f/ds^\. In AQC the interpolation 
from i7 a d(0) to if a d(l) can be chosen at will, in particular 
so as to satisfy the above conditions on H al \ . This shows 
that closed-system AQC is resilient against control errors 
which cause H a ^{s) to deviate from its intended path, as 
long as these do not modify the end point i/ a d(l). This is 
a form of inherent fault tolerance to control errors which 



is not shared by the circuit model 17 1. 

Open system evolution. — A description in terms of _ff a d 
alone neglects the fact that in reality the adiabatic quan- 
tum computer system is never perfectly isolated. The ac- 
tual Hamiltonian is H(t) — Hs(t)^)Is +Is<&Hb +Hsb, 
where I denotes the identity operator, Hs = H a< \ + He 
(Hb) acts on the system (bath) alone, Hc(t) is a control 
Hamiltonian, and Hsb = J^a^a <S> B a , where S a (B a ) 
acts on the system (bath). The role of Hq is to imple- 
ment a DD procedure. The total propagator is U(t) = 



2 



Texp[— i J Q H[t')dt'], where T denotes time ordering. 
The time evolved system state is ps(t) = Tr B p[t), where 
p[t) = U(t)p(0)U(ty is the joint system-bath state. Be- 
low I explain how to choose Hc[t) so that 



[H ad [t),H c [t')]=0 Vi.t'. 



(2) 



It is this condition that will allow application of DD with- 
out interfering with the adiabatic evolution. Consider 
the uncoupled setting Hgg = 0, to be denoted by the su- 
perscript 0. The ideal, noise-free adiabatic system state 
is Pg ad (i) = |0ad(t))(0ad(£)|- Because the adiabatic, 
control, and bath Hamiltonians all commute we have 
p°[t) = P°s(t)®Pc(t)®PB(t), where p° s [t) = |V(i)> W)l 
[Pc(t) = \^c(t))(ipc(t)\] is the actual system evolution 
under H^d [He], and p%[t) is the bath state evolved un- 
der H B . Let p^t) = p° SM (t) ® p% (t) ® p° (t) denote 
the "ideal adiabatic joint state," with purely adiabatic 
evolution of the first factor. Note that /95(C)) = Ps a( j(0)- 
General error bound. — Let d (5) denote distances in 
the joint (system) Hilbert space. To quantify the devia- 
tion of the actual evolution from the desired one, let: 

5 S = D[ps(T),p Sad (T)), d u = D[p[T),p°[T)} 
4d = D[p° (T) , p® d (T)} = <5 ad , d tot = D[p(T), p° d (T)}. 

The overall objective is to minimize the distance 5s be- 
tween the actual system state and the ideal, noise-free 
adiabatic system state. The distance between the un- 
coupled joint state and the ideal adiabatic joint state is 
d a d, which equals 5 ao \ since ||Acg)-B||i = ||A||i||B||i and 
IIPbIIi = II Pell 1 = 1- The "decoupling distance" is du- 
the distance between the joint state in the coupled and 
uncoupled settings. Minimization of this distance is the 
target of the DD procedure. Finally, d to t is the distance 
between the actual and ideal joint states. 

Because taking the partial trace can only decrease the 
distance between states 0] , we have 5s < dtot ■ Using the 
triangle inequality we have dtot < du + d a d- Therefore: 



Ss < du + 



ad • 



(3) 



This key inequality shows that the total system error 
is bounded above by the sum of two errors: (i) due to 
the system-bath interaction in the presence of decou- 
pling (cZd); (ii) due to the deviations from adiabaticity 
in the closed system (d a d)- I shall present a procedure 
intended to minimize dr> jointly with d a d- This is an 
optimization problem: generically decoherence (closed- 
system adiabaticity) worsens (improves) with increasing 
T. 

Dynamical decoupling. — I now show how to minimize 
the decoupling error To do so I propose to apply 
strong and fast dynamical decoupling (DD) pulses to the 
system on top of the adiabatic evolution. It is conve- 
nient to first transform to an interaction picture defined 
by Had + H B , i.e., U[t) = f/ad(f) ® U B [t)U[t), where 



U x (t) = Tex.p[-if*H x (t')dt'], X e {ad,B}. Then U 
satisfies the Schrodinger equation dU/dt = —iHU, with 
H = UI® Ul d [H G + HsbPb ® £4d = H c + H SB , where 
the second equality required Eq. (|2|). Define an effec- 
tive "error Hamiltonian" H c g(t) via U(t) — e - ** 5 "*^'), 
which can be conveniently evaluated using the Mag- 
nus expansion [19] . Now consider a sequence of non- 
overlapping control Hamiltonians ff^ [t) applied for du- 
ration w (pulse width) at pulse intervals r, i.e., Hc[t) = 
for tfc < i < tk+i — w and Hc(t) = H^ for tk+i — w < 
t < tk+i, where t k — k[r + w), k <E TLk- The se- 
quence 

{^ddI^o 1 defines a "DD protocol" with cycle 
time T c — K(t + w) and unitary pulses P k generated 



by H[t) = ffW + H SB , t k+1 



the "ideal pulse limit" w 
pling group" Q = {G k = 



w < t < tk+i- In 
one defines the "decou- 



P K -i---Pk+iPk} 



K-l 

k=0 



such 



that Go = Is- Then the total propagator becomes 
U(T C ) = Y\k=^M-^(G\Hs B G k )] = e- iT ° H i*, where 
H^g denotes the resulting effective Hamiltonian, with 

Magnus series H^ = J2T= [8]. To lowest order: 



K-l 



H. 



id(0) 



cir 



SB G k =Ug(H SB ). (4) 



k=0 



In the limit r — > one has = H^°\ so that by 
properly choosing Q one can effectively eliminate Hsb- 

Returning to non-ideal [w > 0) pulses, we have shown 
by use of ||L4,B]||i < 2||A||||5||i and the Dyson ex- 
pansion that minimization of the "error phase" $(T) = 
T\\H c g (T)l|implies minimization of the decoupling dis- 
tance du [20(: 



d D < min[l,(e* - l)/2] 
< $ if $ < 1. 



(5) 



For single-qubit systems we and others have shown that 
concatenated DD pulse sequences can decrease $ expo- 
nentially in the number of concatenation levels [21j . Here 
I focus on periodic pulse sequences for simplicity. In pe- 
riodic DD (PDD) one repeatedly applies the DD protocol 

{^ddI^To 1 t0 the system, i.e., H G (t + IK) = H c (t), I G 
Z/,. The total time is thus T = L(t + w), where the total 
number of pulses is L and the number of cycles is L/K. A 
calculation of the total error phase $(T) proceeds in two 
steps. First we find an upper bound 0/ on <&;(T C ) for the 
Ith cycle, using the Magnus expansion. Then we upper 
bound $(T) by [L/K) max, 6,. Let J = \\H S b\\ (system- 
bath coupling strength), (3 = \\H &A + H B \\ < fis + Pb, 
where /3s = ||-ff a d|| and (3 B — \\H B \\, and a = 0[1) a 
constant. A worst case analysis yields (l8| : 

, . a(JT) 2 JTw ,exp(2(3T c ) - 1 N , , 



3 



This bound is valid as long the third term is < JT and the 
Magnus series is absolutely convergent over each cycle, a 
sufficient condition for which is JT C 

Joint AQC-DD optimization. — Recall Eq. ([T]) for 
closed system adiabaticity. The given and fixed parame- 
ters of the problem are J, Ao, and z (or £). The task is 
to ensure that each of the terms in Eq. ([6]) vanishes as a 
function of n. I show in 22] that if r and w scale as 



/A , 



-(2C+ei+e 2 ) 



M (7) 



with ei > 1 and €2 > 0, then 



du < ( J/A ) 2 n- £1 + n-^ + (J/A > 



(8) 



which is arbitrarily small in the large n limit. Combining 
this with the bounds above (5 a( j < r -2 or 5 a <j < r~ N ) 
and inequality (J3j) , it follows that for an AQC algorithm 
with time scaling as T = L(t + w) ~ Aq the total 
error 8s can be made arbitrarily small. This is the first 
main result of this work: using PDD with properly chosen 
parameters we can obtain arbitrarily accurate AQC. 

However, there is a shortcoming: the pulse intervals 
and widths must shrink with n as a power law, with an 
exponent dictated by the dynamical critical exponent z 
of the model [Eq. 0]. I expect that this can be remedied 
by employing concatenated DD [IH, [2l[ . 

Seamless AQC-DD. — The entire analysis relies so far 
on the "non-interference" condition ((2|). When can it 
be satisfied? Fortunately, the general background theory 
was worked out in [l(| , though without any reference 
to AQC. I review this theory and make the connection 
to AQC explicit. The decoupling group Q induces a de- 
composition of the system Hilbcrt space TLs via its group 
algebra CQ and its commutant CQ', as follows: 

H s ~0C" J ®C dj , (9) 
,7 

eg s* i nj ® M dj , eg' s* M nj ® i dj . (10) 

Here nj and dj are, respectively, the multiplicity and 
dimension of the Jth irreducible representation (irrep) 
of the unitary representation chosen for £7, while In and 
Mn are, respectively, the N x N identity matrix and un- 
specified complex- valued N x N matrices. The adiabatic 
state is encoded into (one of) the left factors Cj = C™ ' , 
i.e., each such factor (with J fixed) represents an rij- 
dimensional code Cj storing log d nj qudits. The DD 
pulses act on the right factors. As shown in the dy- 
namically decoupled evolution on each factor (code) Cj 
will be noiseless in the ideal limit w, r — ► iff Hg(S a ) = 
(B j ^J,alnj <8> Idj for all system operators S a in Hsb, 



id(0) 



©J [{In., ® hj)] S ® E a 



whence H t 

Thus, assuming the latter condition is met, under the 
action of DD the action of flfjj °^ on the code Cj is pro- 
portional to I n , 7 i-e., is harmless. Quantum logic, or 



AQC, is enacted by the elements of Cg' . Dynamical de- 
coupling operations are enacted via the elements of CQ. 
Condition ([2]) is satisfied because [Cg,Cg'\ = 0. 

Stabilizer decoupling. — An important example of the 
general C(//C(/'construction is when g is the stabilizer 
of a quantum error correcting code and the commutant 
is the normalizer Af of the code 



12] 



Because a stabilizer 
group is Abelian its irreps are all one-dimensional. A 
stabilizer code encoding n qubits into rij = k has n — k 
generators, each of which has eigenvalues ±1. Then J 
runs over the 2"~ fc different binary vectors of eigenval- 
ues, meaning that Hs — © j={±i ±1} C- 2 , and each of 
the subspaces in the sum is a valid code Cj. Here the el- 
ements of Af are viewed as Hamiltonians. For this reason 
only the encoded single-qubit normalizer operations are 
required; encoded two-body interactions are constructed 
as tensor products of single-qubit ones. 

Energy- gap protection. — Application of DD pulses is 
the main mechanism I propose for protection of AQC, 
but it has a shortcoming as noted above. Fortunately, 
the formulation presented here easily accommodates the 
AQC energy-gap protection strategy proposed in [l3j ]. 
which can be viewed as adding another layer of pro- 
tection for dealing with finite-resource-DD. Namely, if 
the decoupling group g is also a stabilizer group for 
code Cj, then for each Pauli error S a in Hsb there is 

0, 

can then add 

e eg 



at least one element Pj € g such that {Pj,S a } 



[IJ. We 



and otherwise [Pj^S a ] 

an energy penalty term Hp = — -EpXwJi 



to Hs, where Ep > 

coupling means that i^g 

o-id(i) 



o(l) 



is the penalty. Imperfect de- 
7^ 0. To lowest order, 
and an "erred state" will be 

M 



P, 



of the form = S a \ip), where 

Vj. Then H P \1>£) = {[« - (K- 1)] (K - 1)E P } |^), 
where a is the number of stabilizer elements that anti- 
commutc with S a . Thus \tp^) is an eigenstate of Hp 
and has a(K — l)Ep more energy than any state in the 
code space. Ref. hj showed, using a Markovian model 
of qubits coupled to a photon bath, the important re- 
sult that this energy gap for erred states implies that 
the temperature need only shrink logarithmically rather 
than polynomially in the problem size. However, note 
that to deal with generic system-bath interactions both 
the stabilizer and normalizer elements must involve k- 
local interactions, with k > 2 [l3j ], 

2-local decoherence-resistant universal AQC. — First 
recall a recent universality result. The following sim- 
ple 2-local Hamiltonian allows for universal AQC [23| : 

With this all the tools have been assembled to demon- 
strate the second main result of this work: a stabi- 
lizer decoupling procedure against 1-local noise that uses 
only 2-local interactions. By 1-local noise I mean the 
main nemesis of quantum computing, namely the lin- 



ear decoherence model: Hgg 



4 



where {B"} are arbitrary bath operators. To beat 
Hgg, use the Abelian "universal decoupling group" [8| 
G nni = {I,X,Y,Z}, where X(Y,Z) = ®'; 
It is simple to verify that ^lg lmi (Hgg) = 0. As 
noted in Ref. Q, un i is the stabilizer of an [[n, n — 
2,2]] stabilizer code C, whose codewords are {\ip x ) = 
(\x) + | not a;)) /V2}, where x is an even- weight binary 
string of length n, with n even. For example, for 
u = 4we find: |00) L = (|0000> + |1111))/V2, |10) L = 
(|0011) + |1100)) A/2, |01) L = (|0101) + |1010)) /V2, 
\11) L = (|1001) + |0110)) /V2. Now universal AQC over 
C can be implemented using 2- local Hamiltonians. To 
compute over C we replace each Pauli matrix in H^ lv 
by its encoded partner. Encoded single-qubit operations 
for C are the 2-local Xj = crf(Tj +1 and Zj — (Jj +1 cr^, 
where j = l,...,n— 2. The 2-local interactions crfuj 
and of cj appearing in _ff a( j are replaced by the 2-local 
XiXj = crf +1 crj +1 and ZiZj = of +1 cJ +1 - Thus we see 
that universal AQC can be combined with DD using only 
2-local afa^ and afa* interactions over C. 

Examples of promising QC implementations where X, 
Z (as pulses for DD) and ofcrj, of crj (as Hamiltonians for 
AQC) are available and controllable, are systems includ- 
ing capacitive coupling of flux qubits |24| and spin models 
implemented with polar molecules [25| . Also note that in 
principle, as discussed above, we can create an additional 
energy gap [13l ] against single-qubit errors by adding a 
penalty term Hp = —Ep(X + Y + Z) to the system 
Hamiltonian. However, Hp is an ro-local interaction. 

Conclusions and outlook. — Using a combination of 
various tools in the arsenal of decoherence control I have 
shown how to protect AQC against decoherence. While 
I believe that the methods proposed here should signifi- 
cantly contribute towards the viability and robustness of 
AQC, what is still missing is a threshold theorem for fault 
tolerant AQC. This will most likely require the incorpo- 
ration of feedback, in order to correct DD pulse imper- 
fections and other control noise [l3]- One possibility for 
doing so might be to perform syndrome measurements on 
the commutant factor [C dj in Eq. ©] as in recent circuit- 
model fault tolerance work using subsystems codes [7|. 
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